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Here we report from our theoretical studies that in biased bilayer graphene, one can induce phase 
transitions from an incompressible fractional quantum Hall state to a compressible state by tuning 
the bandgap at a given electron density. The nature of such phase transitions is different for weak 
and strong inter-layer coupling. Although for strong coupling more levels interact there are lesser 
number of transitions than for the weak coupling case. The intriguing scenario of tunable phase 
transitions in the fractional quantum Hall states is unique to bilayer graphene and never before 
existed in conventional semiconductor systems. 



The unconventional quantum Hall effect in monolayer 
graphene, whose experimental observation [2| unleashed 
quite unprecedented interest in this system [3(, reflects 
the unique behavior of massless Dirac fermions in a mag- 
netic field [J, Q . In bilayer graphene this effect confirms 
the presence of massive chiral quasiparticles [6| . An im- 
portant characteristic of bilayer graphene is that it is a 
semiconductor with a tunable bandgap between the va- 
lence and conduction bands 7]. This modifies the Lan- 
dau level spectrum and influences the role of long-range 
Coulomb interactions Q. Here we report that the frac- 
tional quantum Hall effect (FQHE), a distinct signature 
of interacting electrons in the system [9|, LL(| is very sensi- 
tive to the interlayer coupling strength and the bias volt- 
age. We propose that by tuning the bias voltage one can 
induce phase transitions from an incompressible state to 
a compressible state at a given gate voltage. The bilayer 
graphene system shows quite different properties for weak 
and strong inter-layer coupling. For a weak coupling 
the energy spectrum as a function of bias voltage shows 
a set of anti-crossings, resulting in transitions from the 
FQHE state to a compressible state. At strong coupling 
there is a strong interaction between many energy lev- 
els, which finally results in only a few phase transitions. 
This interesting scenario of tunable phase transitions in 
the FQH states is unique to bilayer graphene. In con- 
ventional semiconductor systems the type of phase tran- 
sitions discussed below was never reported. The FQHE 
in monolayer graphene was in fact, studied theoretically 
by us [ll[ and subsequent experiments confirmed the ex- 
istence of that effect in suspended monolayer graphene 
samples [12|. No such studies have been reported on bi- 
layer graphene. 

We assume that the bilayer graphene consists of two 
coupled graphene layers with the Bernal stacking ar- 
rangement. Our main concern then is the coupling be- 
tween atoms of sublattice A of the lower layer and atoms 
of sublattice B' of the upper layer. The single-particle 
levels have two-fold spin degeneracy and two-fold valley 
degeneracy, which can be lifted in the many-particle sys- 



tems at relatively large magnetic fields [13|. The valley 
degeneracy is also lifted under an applied bias voltage 
6] . Considering only one spin direction, we describe the 
state of the system in terms of the four-component spinor 
{i\) A ,vb B ,vbB'-,^A') T for valley K and {^) B ','^A','>P a ,^b) T 
for valley K' . Here subindices A, B and A', B' correspond 
to lower and upper layers respectively. The strength of 
inter-layer coupling is described in terms of the inter-layer 
hopping integral, t. In a biased bilayer graphene the bias 
potential is introduced as the potential difference, AC/, 
between the upper and lower layers. The Hamiltonian 
of the biased bilayer system in a perpendicular magnetic 
field then takes the form [6| 
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— n x s- m y , vi = p + eA/c, p is the two- 
dimensional electron momentum, A is the vector poten- 
tial, v F sw 10 6 m/s is the fermi velocity, and £ = + (K 
valley) or — (K' valley). 

In a perpendicular magnetic field the Hamiltonian 
(fTj) generates a discrete Landau level energy spectrum. 
The corresponding eigenfunctions can be expressed in 
terms of the conventional nonrelativistic Landau func- 
tions. The electron states in sublattices A and A' are 
written in terms of the ro-th Landau functions, while the 
electron states in sublattices B and B' are described by 
the \n — 1| and n + 1 Landau functions, respectively. 
Therefore the Landau states in bilayer graphene can be 
described as a mixture of n, n + 1, and n — 1 nonrela- 
tivistic Landau functions belonging to different sublat- 
tices [7j . This mixture, for a given value of n, results in 
four different Landau levels. The Landau level energies, 
e corresponding to the index n can be found from the 
following equation [7] 

\e + tff - 2(n + l)][(e - £<5) 2 - 2n] = (e 2 - S 2 )t 2 , 

(2) 



where S = AU/2 and all energies are expressed in units 
of hv F /£ . Here £ Q = {h/eB)? is the magnetic length. 

We now introduce a labeling scheme for Landau levels 
in bilayer graphene. From Eq. @, we see that for each 
value of n (= 0,1,2,...) and in each valley there are 
four solutions, i.e., four Landau levels. Usually, the two 
lower Landau levels have negative energies and the two 
upper Landau levels have positive energies. Then each of 
the four Landau levels can be labeled as n\ , where i = 
—2, —1, 1, 2 is the label of the Landau level corresponding 
to the solution of Eq. @ for a given value of n in the 
ascending order. 

For a partially occupied Landau level the properties 
of the system, e.g., the ground state and excitations, are 
completely determined by the inter-electron interactions, 
which can be expressed by Haldane's pseudopotentials, 
V m , [3 (energies of two electrons with relative angu- 
lar momentum to). In a graphene bilayer the Haldane 
pseudopotentials in a Landau level with index n and the 
energy e have the form 



yip.) 
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where L m {x) are the Laguerre polynomials, V(q) = 
2ire 2 / (n£ q) is the Coulomb interaction in the momen- 
tum space, k is the dielectric constant, and F n ^ s (q) are 
the corresponding form factors 
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The form factors of bilayer graphene [Eq. ([JJ] are 
clearly different from the corresponding ones for a mono- 
layer graphene. In the latter case, the form factor of the 
n = Landau level is the same as that of conventional 
non-relativistic electrons [ll|, LLIIj Fo(q) = Lq. The form 
factors of higher Landau levels are determined by the 
mixture of L n and L n+ \ terms. In bilayer graphene the 
form factors of the n = Landau level are mixtures of 
the Lq and L\ terms and are different from that in the 
non-relativistic case. There is one special Landau level in 
bilayer graphene with index n — 0, whose properties are 
completely identical to that of the non-relativistic n = 
Landau level. It is clear from Eq. © that for n = 
there is a Landau level with energy e = £<5. This energy 
does not depend on the coupling between the layers, t. 
The form factor of this Landau level is exactly equal to 
the form factor of a non-relativistic system of the n = 
Landau level, F n= o, e= ^s = Lq. Therefore, all many-body 
properties of a bilayer system in the n = 0, e = £6 Landau 
level are completely identical to those of a non-relativistic 
conventional system in the n = Landau level. 



For Landau levels with higher indices, the form fac- 
tor is a mixture of three different functions, L n , L n —±, 
and Lji+i- Therefore, in general, the strength of inter- 
electron interactions in bilayer graphene is strongly mod- 
ified as compared to its value in monolayer graphene. To 
address the effects of these modifications on the prop- 
erties of the many-electron system in bilayer graphene 
we investigate fractional filling factors corresponding to 
the FQHE [9j. We treat the many-electron system at 
various fractional fillin g fa ctors numerically within the 
spherical geometry [ll], [lj| . The radius of our spehere is 
R = \S£ , where 2S is the number of magnetic fluxes 
through the sphere in units of the flux quanta. The 
single-electron states are characterized by the angular 
momentum, 5, and its z component, S z . For the many- 
electron system the corresponding states are classified 
by the total angular momentum L and its z component, 
while the energy of the state depends only on L [lj| . A 
given fractional filling of the Landau level is determined 
by a special relation between the number of electrons N 
and the radius of the sphere R. For example, the |- 
FQHE state is realized at S = (|)(JV - 1), while the 
|-FQHE state corresponds to the relation S = {\)N — 2. 
With the Haldane pseudopotentials [Eq. (J3H we deter- 
mine the interaction Hamiltonian matrix [16j | and then 
calculate a few lowest eigenvalues and eigenvectors of this 
matrix. The FQHE states are obtained when the ground 
state of the system is an incompressible liquid, the energy 
spectrum of which has a finite many-body gap [9|, [l0( . 



We begin with the celebrated ^-FQHE [lOj, corre- 
sponding to the filling factor v — \. The behavior of 
the Landau level spectra as a function of the bias voltage 
and for different values of t are displayed in Fig. Q] (only 
the Landau levels with positive energies are shown). A 
similar behavior is valid for other FQHE filling factors, 
e.g., for v = g [r7| . Figure [1] clearly illustrates that the 
FQHE can be observed in all n = Landau levels with 
the strongest FQHE being in the second n = Landau 
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levels of both valleys, i.e., O2 and 

We found an interesting behavior in the n = 1 Landau 
levels. There are four such levels with positive energy; 
two per each valley. The FQHE in these levels shows dif- 
ferent properties depending on the strength of t. For all 
parameters of the system there is no FQHE in the Landau 
level l'fK At small values of t, t < 150 meV, (Fig. Ufa)) 
the system clearly shows few anti-clossings accompanied 
by the transitions from the FQHE incompressible state 
to a state without the FQHE. There is one such transi- 
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tion for levels 1 2 and 1\ , but there are two transitions 

, (+) 



in level 1\ , corresponding to two anti-crossings in this 

level. Thus, in level l^ +) and small AU, the FQHE is 
present but disappears at larger values of AU. It reap- 
pears at very large values of Af (~ 400 meV). With in- 
creasing t the two anti-crossings in level 1^ merge (see 
Fig.[TJb)) and finally disappear (Fig.QJc)). At large val- 
ues of t, t > 150 meV, there are only two anti-crossings 
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FIG. 1: A few lowest Landau levels of the conduction band 
(for two valleys) as a function of the bias potential, AU, for 
different values of inter-layer coupling: (a) t = 30 meV (b) t = 
150 meV and (c) t — 300 meV and a magnetic field of 15 Tesla. 
The numbers next to the curves denote the corresponding 
Landau levels. The Landau levels where the i-FQHE can be 
observed are drawn as blue and green filled dots. The green 
dots correspond to the Landau levels where the FQHE states 
are identical to that of a monolayer of graphene. The red dots 
represent Landau levels with weak |-FQHE and the open dots 
for those where the FQHE is absent. In (a), the dashed lines 
labeled by (i) and (ii) illustrate two situations: (i) under a 
constant gate voltage and variable bias potential; (ii) under a 
constant bias potential and variable gate voltage. 
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(Fig.QJc)) in 1\~' and 1 2 Landau levels. At such large 
values of t, the anti-crossings cannot be considered as in- 
teraction between two 'crossing' levels, but as a result of 
strong interaction between all (four) levels of the two 
layers of bilayer graphene. It is important that such 
strong interaction between the levels does not destroy 
the FQHE, but shows well-defined regions with strong 
FQHE. For weak coupling between graphene monolay- 
ers, i.e. for a small t, transitions from the |-FQHE state 
to a non-FQHE state can be understood in terms of the 
anticrossing of n = 1 and n = 2 Landau levels of the 
monolayers. For monolayers, the FQHE can be observed 
only in the n = and 1 Landau levels but not for n = 2 



ll[ . The levels without the FQHE in Fig. Q] (c) then cor- 
respond to n — 2 of one of the monolayers. For large t, 
i.e. for strong coupling, such a simple description is how- 
ever inadequate. The properties of the n = 1 levels have 
important implications for possible experimental obser- 
vations of this unique behavior (Fig.[TJa)): 

(i) By applying a gate voltage the electron density can 
be tuned so that the first four Landau levels are com- 
pletely occupied and the next Landau level is partially oc- 
cupied with the FQHE filling factor, for example, v = |. 
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FIG. 2: Low-energy excitation spectra of the i-FQHE states 
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(eight electrons) in the Landau levels (a) lL. , and (b) 1^. , 
shown for different values of the bias potential (the numbers 
next to the lines are the values of AU in meV). The system is 
fully spin-polarized. The inter-layer hopping integral is set to 
30 meV and the magnetic field is 15 Tesla. The flux quanta 
is 2S = 21. The solid dot at L — depicts the ground state. 



The energy unit is e c 
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Following Fig. Ufa), this means that the 0^, o[ +) , 
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Og Landau levels are fully occupied, while the lj Lan- 
dau level has a filling factor |. Then, by varying AU from 
a small value, e.g., 10 meV, to a larger value, e.g., 200 
meV, one can observe the disappearance of the FQHE 
(line (i) inFig.QIa)). 

(ii) The bias voltage is kept fixed at a large value, 
e.g., AU = 300 meV. Then by varying the gate voltage 
and thus increasing the electron density, one can observe 
the disappearance and reappearance of the 1/3-FQHE 
in higher Landau levels (line (ii) in Fig. [TJa)), when the 
filling factors of the corresponding Landau levels are 1/3. 

The collapse of the FQHE gap corresponding to the 
appearence of anticrossing of the n — 1 Landau levels, 
is illustrated in Fig. [5] The FQHE gap has a monotonic 
dependence on the bias voltage. In the anticrossing re- 
gion the gap disappears for the lower n = 1 Landau level 
(Fig. |5Ji) and reappears for the higher n = 1 Landau 
level (Fig. |2b)- The evolution of the energy spectra of 
the incompressible liquid is found to be similar for other 
filling factors (such as v = | [17|]). This behavior was 
never before observed in the FQHE of conventional two- 
dimensional electron systems. 

The strength of the FQHE, i.e., the magnitude of the 
excitation gap, depends on the bias voltage and the inter- 
layer hopping integral. In Fig. [3l this dependence is 
shown for i-FQHE in different Landau levels as a func- 
tion of t. In accordance with the properties of Haldane 
pseudopotentials, the excitation gap of the 0^ Landau 
levels does not depend on the bias voltage and on the 
inter-layer hopping integral. The corresponding gap re- 
mains constant and is equal to the gap of the FQHE in 
a single layer of graphene in the n — Landau level. For 
t = the two layers of graphene are decoupled and the 
bilayer system becomes identical to a monolayer with ad- 
ditional double degeneracy. This property is clearly seen 
in Fig. |3l where for t = there are only two doubly de- 
generate FQHE gaps, corresponding to n = and n = 1 
single layer Landau levels. 

At the zero bias voltage the system has two-fold valley 
degeneracy, which is lifted at finite values of AU. At a 
small bias voltage, AU = 10 meV, the levels belonging 
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FIG. 3: The FQHE gaps are shown for different Landau levels. 
The labels next to the lines correspond to the labeling of 
Landau levels shown in Fig. [l] v = i-FQHE (eight electron) 
for (a) AU = 10 meV, and (b) AU = 300 meV. All systems 
are fully spin polarized. 



to two valleys are almost degenerate, which results in al- 
most the same FQHE gaps of the corresponding levels. 
At the same time the FQHE gaps of 0[ and 0^ levels, 
which are degenerate at the zero bias voltage, are differ- 
ent. The origin of these levels is the following: At the 
zero bias voltage there is a four-fold degenerate Landau 
level with zero energy (0 1 +) , ()£ -) , 0^, and 0^). At a 
finite bias voltage, two of the levels have positive energies 
(shown in Fig. Q]) and the other two levels have negative 
energies. At small values of AU, the wavefunctions cor- 
responding to the levels 0^ and 0_^ of valley K have 
the form (0, 0, 0, (f> ) and (^ , 0, 0, {t/2^)^), respectively. 
Here 4> n are n-th 'nonrelativistic' Landau functions and 
t is in units of hv F /£ . The corresponding form factors 

F(q), are L for level o[ +) and L °\f t if^ for level 0^ } . 
Although the energies of these levels are almost the same 
the form factors and hence the gaps are quite different. 
In Fig. [3ja) this difference is clearly visible. The depen- 
dence of the gap of FQHE of level 0\~' on parameter t is 



nonmonotonic. At t = the form factor of level 0^ is 
L and the FQHE gap is exactly equal to the FQHE gap 
in the n = Landau level of a single graphene layer. At 
t = 2?hv F /£ the form factor is L °+ £l and the FQHE 
gap is equal to that in the n — 1 Landau level of a single 
graphene layer. This point corresponds to the maximum 
in Fig. EJa). At a large bias voltage, ALT = 300 meV 
(Fig. [3}d), the FQHE gaps show mainly monotonic de- 
pendence on the hopping integral. The FQHE gaps of 
the levels 0\ and Q\ , which were quite different in 
Fig. El^a), are now close. There is also disappearance of 
the FQHE in some n = 1 Landau levels, which corre- 
sponds to the anti-crossing behavior in Fig. [TJ Similar 
results are also found for the filling factor v = f (i.e., 
FQHE o non-FQHE transitions occur at the same Lan- 
dau levels and at similar values of AU) |17| . 

To summarize, we have clearly demonstrated that bi- 
layer graphene in a strong perpendicular magnetic field 
reveals some unique properties, which could allow novel 
transitions from the FQHE state to a vanishing FQHE 
state. These transitions occur within the same Landau 
level by varying the bias voltage, i.e, the potential dif- 
ference between the layers. Similarly, we have shown 
that our work on bilayer graphene also results in new 
physics: The transitions FQHE O zero-FQHE, which for 
weak inter-layer coupling can be explained as the result 
of anti-crossing of two levels, also persists in the limit 
of strong coupling, where all levels are strongly coupled. 
We have established here that there is a fundamental 
difference between the two regimes of weak and strong 
coupling in bilayer graphene. The boundary between the 
two regions is determined by the dimensionless parame- 
ter [t/(hv F /£ Q )} ~ 1.5. 
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